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Abstract 

Charged-current cross sections are calculated for quasielastic neutrino 
and antineutrino scattering using a relativistic meson-nucleon model. We 
examine how nuclear-structure effects, such as relativistic random-phase- 
approximation (RPA) corrections and momentum-dependent nucleon self- 
energies, influence the extraction of the axial form factor of the nucleon. RPA 
corrections are important only at low-momentum transfers. In contrast, the 
momentum dependence of the relativistic self-energies changes appreciably 
the value of the axial- mass parameter, Ma, extracted from dipole fits to the 
axial form factor. Using Brookhaven's experimental neutrino spectrum we 
estimate the sensitivity of M^i to various relativistic nuclear-structure effects. 
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I. INTRODUCTION 



The strange-quark content of the nucleon has received considerable attention as a result 
of the measurement of the spin- dependent structure function of the proton by the European 
Muon Collaboration (EMC) 0. Some analyses of the experiment suggest that a large 
portion of the spin of the proton is carried by strange quarks. One can attempt to resolve 
this "spin problem" by studying the strange-quark contribution to the vector (both electric 
and magnetic) and axial-vector form factors of the nucleon. Most likely, it will take a large 
number of measurements to determine all of these form factors separately. Moreover, there 
are important complications from radiative corrections 0, which hinder the extraction of 
strange-quark matrix elements from parity violating electron scattering. Therefore, one 
anticipates a program of several electron experiments || which, combined with neutrino 
scattering data Q, will offer the most accurate strange-quark information. 

Neutral-current neutrino scattering is sensitive to the strange-quark matrix elements of 
the nucleon — especially to the isoscalar component |J. Complications arise, however, 
from the fact that most neutrino experiments measure a combination of elastic scattering 
from free protons plus quasielastic scattering from nucleons bound in nuclei. In the present 
work we examine how nuclear-structure corrections affect the extraction of strange-quark 
information. 

The isoscalar part of the axial-vector form factor of the nucleon is characterized by 
the parameter g s A — the value of the isoscalar strange form factor at zero four-momentum 
transfer. Unfortunately, the extraction of g s A from the Brookhaven National Laboratory 
(BNL) experiment is complicated by low statistics @]. Moreover, there is a strong correlation 
between the extracted value of g s A and the axial mass Ma obtained from dipole fits to the 
axial- vector form factor ||[7]]. Indeed, the world-average value of Ma (1.032 ± 0.036 GeV) 
that the BNL group has used is significantly different from the one [1.09 ± 0.03 (stat)±0.02 
(syst) GeV] extracted later from a charged-current experiment ||. Thus, this large difference 
in Ma is sufficient to change the value of g A extracted from the BNL experiment and the 
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conclusion of a nonzero strange-quark content from our previous neutral-current study 
Therefore determining the precise value of Ma is important for strangeness studies. 

At relatively high momentum transfer, the response of the nuclear target seems to be 
adequately described in a relativistic Fermi gas (RFG) model. Indeed, at a qualitative level, 
a RFG calculation of quasielastic (e, e') longitudinal and transverse responses agrees well 
with finite-nucleus results. However, is the RFG model accurate enough to determine the 
precise value of Ma ? Currently, the BNL experiment claims an Ma value with an error of 
less than 5 percent. Is there any uncertainty comparable to this small error induced from 
nuclear-structure effects? It is the aim of this work to go beyond the RFG response and 
examine the sensitivity of quasielastic neutrino scattering — particularly Ma — to a variety 
of nuclear-structure corrections. 

One such correction arises from long-range RPA correlations. In a recent paper, Singh 
and Oset used a nonrelativistic RPA formalism to study the nuclear response in quasielastic 
neutrino scattering HI. They found RPA corrections to be large only at low- momentum 



transfers. Since most neutrino experiments (TO] are carried out at medium to high momen 



turn transfers, they concluded that the value of Ma extracted from these experiments is 
reliable. 

In a relativistic description of the nuclear target, additional nuclear-structure corrections 
must be considered. According to quantum hadrodynamics (QHD) [|lj]], the saturation of 
nuclear matter arises from a cancellation of strong scalar (a) and vector (u) mean fields. 
The strong scalar field decreases the nucleon mass while the vector field shifts the four- 
momentum of the nucleon in the medium. These relativistic effects were not addressed in 
the nonrelativistic calculation and may provide interesting corrections to the RFG response. 

In a charged- current reaction, the mean-field ground state can be characterized in terms 
of an effective nucleon mass that is reduced, relative to its free-space value, by the presence 
of the strong scalar field. In turn, the RPA response can be modeled from a (tt + p + g') 
residual isovector interaction. In a mean-field approximation, the effective nucleon mass is 
obtained from solving self-consistently the equations of motion at a given baryon density ||11|| . 



In the case of polarized electron scattering, the parity-violating asymmetry was found to be 



sensitive to the in-medium value of the nucleon mass |T2[ . This is an interesting result that 
should be incorporated in neutrino-scattering studies. 

This paper is organized as follows. Section [IT] presents the formalism for the cross section 
of inclusive neutrino scattering in a relativistic impulse approximation and in RPA. Results 



for the charged current cross section are presented in section |HI[ while section [TV| is a 
summary. 



II. FORMALISM 

In a charged- current process, neutrinos (and antineutrinos) interact with nuclei via the 
exchange of charged weak- vector bosons (W ± ) with the resulting production of charged 
leptons (electrons or muons) in the final state. In an inclusive process, where only the 
final leptons are detected, the most general expression for the cross section can be given in 
terms of a time-ordered product of current operators. From this general expression several 
approximations can be made depending on how one treats the ground state of the nucleus 
and its response to the external probe. In Sec. [11 A] , we derive a general formalism for 



the inclusive process and discuss various approximations for the response in a mean-field 



approximation to the ground state. In Sec. |IIB| we discuss the form of the nuclear current 



adopted in the calculation while Sec [II C| contains a detailed description of the relativistic 
random phase approximation. Finally, in Sec. [II D| we discuss how the mean-field ground 
state is modified by the introduction of phenomenological momentum-dependent corrections 
to the nucleon self-energies. 



A. GENERAL FORMALISM 

The scattering process we consider is shown in Fig. [1]. An incoming neutrino with 
momentum k scatters off the nucleus via the exchange of weak-vector bosons producing a 
charged lepton with momentum k! in the final state. The initial and final states of the nucleus 
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are denoted by \if>i(pi)) and \i>f{pf)), respectively. In Born approximation the inclusive cross 
section becomes proportional to the contraction of a leptonic and a hadronic tensor: 

da oc L^W^ (1) 

where the corresponding leptonic (L^) and hadronic (W tJ,v ) tensors are given by 



h' h — h ■ h 1 n 4- h' h zsz ip a ^ u h la h^ 
ft, hj V ft, ft, y^ u T n v n^ -\- It ft, ft, 



(2) 



W" = £ (2vr) 4 5 (4) fe + q-Pf) (^|>(0)|^/) (^/|>(0)|^) ■ (3) 
/ 

Here J M is the weak charge- changing current operator of the nucleus, q = k — k' is the 
momentum transfer to the nucleus, and the plus (minus) sign in L^ u corresponds to an- 
tineutrino (neutrino) scattering. Note that our convention for the antisymmetric tensor is 

£ 0123 = y 

Now we introduce the current-current correlation function, or polarization tensor, as a 



time-ordered product of nuclear currents Jl7j 



iW»(q) = f Ae^ 1 (Vi|T(>(aO J u (0)M) , (4) 



The hadronic tensor, and therefore the cross section, can be directly related to the polar- 
ization tensor. In particular, it is easy to show that the cross section takes the following 
form 

da oc Im (L^U^) . (5) 

This expression is convenient for the evaluation of the inclusive response of a many-body 
system like the nucleus. In particular, various approximations can be made depending on 
how one calculates the ground state of the nucleus and its linear response to the external 
probe [|T3, 18 1 . 



For the many-body current operator J M we assume a simple one-body form: 

J M (x) = ij(x)T^(x), (6) 



where if)(x) is a nucleon-field operator and T M is the weak-interaction vertex to be discussed 
below [see Eq. fll^D]. Meson-exchange currents represent corrections to this one-body form 
and will be ignored throughout this paper. In a mean-field approximation to the nuclear 
ground state the time-ordered product of currents can be evaluated readily using Wick's 
theorem, i.e., 

lWV = \wf Tr[G{p + q) r " G{p) n ' (7) 

where G(p) is the nucleon propagator that will be evaluated in various approximations. 

The simplest approximation that we employ treats the nuclear ground state as a rela- 
tivistic free Fermi gas. Here the nuclear response consists of the excitation of particle-hole 
pairs subject to the constraints imposed by energy-momentum conservation and the Pauli 
principle. The nucleon propagator differs from the well known Feynman propagator only 
because of a finite-density correction arising from the filled Fermi sea [|TT|1 , 



G°(p) = (tf + M) 



^^ + |>-W-IP|) • (8) 



Note that we have introduced the Fermi momentum kp and the free (on-shell) energy E p = 
y/p' 2 + M 2 . We call this approximation "impulse with M" in order to distinguish it from 
the self-consistent impulse approximation with an effective mass M* that we now address. 

One can improve the free Fermi-gas description by taking into account, at least at the 
mean-field level, the interaction between the nucleons in the nucleus. In a mean-field- 
theory approximation (MFT) to the Walecka model the propagation of a nucleon through 
the medium is modified by the presence of constant scalar and vector mean-fields. These 
potentials induce a shift in the mass and in the energy of a particle in the medium and give 



rise to a self-consistent nucleon propagator |TT 
G*(p) = (f + M 



^ *(p5 - E*)0{h F - |p|) 



(9) 



p*2 _ m*2 + ie ft* 

where the effective mass and energy are shifted from their free-space value by the scalar (S) 
and timelike component of the mean fields, 



M* = M + S ; E% = \/p 2 + M* 2 ; p*^ = (p° - V , p) . (10) 

These changes in the nucleon propagator induce a corresponding change in the polarization 
tensor which is now written, 

•H£f = / (0iTr[G*(p + g)r^*(p)n . (11) 

Note, in computing the response we integrate over the four-momentum of the nucleons, and 
the contribution from the constant vector potential can be eliminated by a simple change 
of variables. [This will not happen once momentum-dependent corrections are incorporated 
into the mean fields (see Sec. [II D| ) .] Formally the mean- field response is identical to that of 
a relativistic Fermi gas of nucleons with an effective mass M*. We refer to this calculation 
as "impulse with M*" . 



B. IMPULSE APPROXIMATION 

We start this section by writing the inclusive cross section — per neutron — for the 
charge- changing process as 



d 2 a G|cos 2 c |k'| 



Im(L^n^). (12) 



dVL^dE^ 327T 3 pE u 

Here p = k F /3ir 2 is the neutron (or proton) density of the system, 6 C the Cabbibo angle 
(cos 2 # c = 0.95), Gp is the Fermi constant, k' the three-momentum of the outgoing lepton, 
and E v the energy of the incoming neutrino (or antineutrino). 

In the impulse approximation the interaction between the incoming neutrino and a target 
nucleon is assumed to be the same as in free space. Hence, we employ a charge-changing 
current operator with single-nucleon form factors parameterized from on-shell data. That is 
(suppressing isospin labels), 

T»(q) = F x {Q*)r + iF 2 {Q 2 )^^ - G A (Q 2 )^ + F p (Q 2 )q^ , (Q 2 = q 2 - q 2 ) (13) 

The form factors Fi,F%, Ga and F p are given in Appendix 1. The pseudoscalar form factor 
F p is constructed from PCAC and its contribution is suppressed by the small lepton mass. 



Since T M has been expressed in terms of vector, tensor, axial-vector, and pseudoscalar 
vertices, the inclusive cross section requires the evaluation of a large set of nuclear-response 
functions. These are conveniently separated in the following way (note that the subscripts 



indicate the vertices involved): 

IC = -i / -0^Tr[G( P + q)rG(p)Y) , (14) 

n£ = -tj ?AG(P + q) G(p) -^-] , (15) 

Kt=~i J ^^(P + ^-fGip)^^.) , (16) 

IC = "» / (Si Tr t G G° + ff) 7* G (p) tV] = -•V" 00 QJl m , (17) 

IC = -i J Tr[G(p + g) 7 V C?(p) 7 V] = IC + , (18) 

n£ = / Tr ( G (p + ?) 7 M 7 5 G(p) (-<f h 5 ] = g"g"n<* , (19) 

°» = _i / Tr[G °° + g) q ^ G{p) { ~ qU)l5] = q " qUUpp • (20) 



The various components of the polarization tensor can be computed in the Fermi-gas limit or 
in the mean-field approximation by using either G°(p) [Eq. @] or G*(p) [Eq. @] respectively. 
The imaginary parts of all these polarizations have been calculated analytically following 
the Ref. []20| and are given in Appendix 2. 

The polarizations IC 11^, and II^ are only sensitive to the Lorentz- vector part of the 
weak current and, thus, satisfy current conservation, g M IF^ = V^ v q v = 0. The conservation 
of the vector current plus Lorentz covariance imply that only two components of each of 
these polarizations are independent. These have been chosen to be the longitudinal and 
transverse components which are defined, for example in the the case of IC as 

n L = n 00 - n 11 - -— n 00 (21) 
n T = 21T 22 = 211 33 . (22) 

vv — vv vv \ J 

Here we have assumed a coordinate system with the x = 1 axis along the direction of 
the three-momentum transfer q. Using Lorentz covariance one can isolate the additional 
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responses that arise from the other components of the current and obtain the following 
invariant amplitude 

-Z^rP 17 = LlRl + LtRt ± L va RvA + LaRa + L p Rp . (23) 

The nuclear-structure information is fully contained in the various response functions which 
have been defined in terms of the above polarization tensors: 

R L = {Ft + G\) Yi L vv + 2F 1 F 2 l£ + Fi n£ , (24) 

Rt = \ [(Ff + G 2 A ) C + 2F,F 2 l£ + Fi IQ , (25) 

R P = 2G A F p U ap + F 2 U pp , (26) 

R A = G A T1 A , (27) 

M* 

Rva = (F 1 + F 2 — ) G A |q| ll va . (28) 

These response functions are multiplied by appropriate kinematical factors that could, in 
principle, be used to separate the individual responses 

q 2 Am 2 qn , on m 2 , 

L L = -- 2 L , «P (AE V - g - , (29) 

L A = 8(q 2 - mj) , (30) 

L T = -^L 00 - ~^(AE V q° - q 2 + mj) - L A , (31) 

L p = Ami K ~ <? 2 ) = ~\ m l L a , (32) 

ir g 2 (^ + ^) + g m; 

L m = -16 — ; £ . (33) 

q 



Finally, we note that the plus (minus) sign in Eq. fl23|) should be used for neutrino (antineu- 
trino) scattering. 

C. RELATIVISTIC RANDOM PHASE APPROXIMATION 

In the present section we improve the simple particle-hole description of the response 
by incorporating many-body RPA correlations. Many-body correlations can be included 
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by considering the residual interaction between the particle and the hole. For the present 
charge- changing reaction only isovector correlations are important. Yet, there might still 
be important effects associated with the isoscalar mean fields (e.g., M*). Indeed, in a 
recent calculation we have shown that the reduction of the effective nucleon mass in the 
medium results in a quenching of the effective NNir coupling which, in turn, is responsible 
for suppressing the predicted enhancement of the spin-longitudinal to spin-transverse ratio, 
in accordance with experiment |L9 . For the residual isovector interaction we employ a 



simple relativistic generalization of the conventional 7r + p + g' interaction PJ21~1. The 
phenomenological Landau-Migdal parameter g' has been included to simulate the effect 
of repulsive short-range correlations. The RPA correction, AYl^p A , to the polarization 
tensor, IP" in Eq. ([12]), is shown diagrammatically in Fig. ^| (a). The RPA corrections 
are calculated from the dressed propagator, Drpa- This includes an infinite sum of the 
lowest-order (uncorrelated) polarization as illustrated in Fig. |2] (b). Note, for the mean- field 
ground state, the nucleon propagators in the lowest-order polarization already include the 
isoscalar dressing due to the mean fields [see Eq. ([5])]. 

The Lagrangian density describing the isovector component of the NN interaction is 
given by 

£ = 9 P ^\ ■ V>P„ + f P ^ v \ ■ ~ ^757^ • d,7z . (34) 

The p meson has a vector (g p ) as well as a tensor (f p ) coupling to the nucleon. The pa- 
rameters of the model are obtained directly from the Bonn potential fit to NN properties 
and are given by g 2 p /^ = 1.64 and f p /g p = 6.1 |22| (note that our value for g p /4:7i is four 
times larger than the one quoted by the Bonn group simply because of our selection of r/2, 
rather than r, as the isospin vertex). With this form for the interaction Lagrangian the 
NNp vertex becomes (combined with isospin matrices) 



■ v - v "- "T^'"" ' 



while the p— meson propagator is given by 
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-Qixv + q^/m: 
q 2 — m? + ie 



= —2 ^2 I 77 ■ ( 36 ) 



Note that the "gauge" piece q^q v /rn 2 p will not contribute to the RPA response because in 
the mean-field approximation the vector-isovector current is conserved. For the NNn vertex 
we have adopted a pseudovector form with /^/47r = 0.075. The pion, rho, and free nucleon 
masses have been fixed at their experimental values. 

The most uncertain parameter in our calculation is the phenomenological Landau-Migdal 
parameter g' . Traditionally, g' is introduced to regularize the large spin-spin component of 
the isovector interaction. In a relativistic formalism we can incorporate short-range correla- 
tion effects by modifying the pion "propagator" in the following way: 



= — y— — -> IV = — 2—— -9 9iw ■ 37 



Without g' the contribution of the pion to the RPA response would be suppressed by current 
conservation and the small leptonic mass. Therefore, it is through the Landau-Migdal 
parameter g' that the pion mainly contributes to the RPA response (note that the term 
proportional to g' has, both, longitudinal as well as transverse components). With this 
choice for the pion propagator the "elementary" NNn vertex becomes 

T NNn = V • (38) 

m n 

With the above Feynman rules in hand we can now construct the medium-modified rho- 
and pion-mediated interactions. Note that the inclusion of g f , which contains transverse 
as well as longitudinal components, is responsible for p — it mixing. In order to properly 
account for this mixing, Dyson's equation for the propagator must be expanded from a 4 x 4 
to an 8 x 8 matrix equation: 

D RPA = Do + D n D fiPA , (39) 



where we have defined the free (diagonal) propagator matrix 

D n - , (40) 



V° v 7 
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in terms of the p [Eq. (|36])] and g'— modified pion propagator [Eq. (p7|)1. We have also 
introduced the mixed p — tt polarization matrix 



f u n N 



(41) 



with individual components given by 



IT = -i [ *E-T: 
pp 2 J (2vr) 4 

P n *p ™ I (2tt) 4 



n 



[IV 



2M 



ml J (2tt) z 



Tr 



2Mg p 

l 5 rG(ph 5 YG(p + q) 



2M 

G(ph 5 YG(p + q) 



(42) 
(43) 
(44) 



(45) 



The RPA correction to the polarization now takes the following form [see Fig. |2| (a)]: 



An 



RPA 



(n„ n,) n RPA 



(46) 



where U p and II,,- characterize the in-medium mixing — due to particle-hole excitations 
of a charged weak-vector boson with a p or tt meson, respectively, and are given by 



V2J (2tt) 4 

U^ = -V2i^ [ -^rTr 
rrwJ 2tt 4 



^G(p)(Y- i 4^ 1 )G(p + q) 



2Mg p 
r"G(p) 7 V G(p + q) 



(47) 
(48) 



An RPA calculation of the inclusive response uses the same expression for the cross section 
as in Eq. ([12]) with the replacement: 



-> n& A = it- + ait^ 



(49) 



In an impulse (or uncorrelated) description of the response the cross section is only 
sensitive to the imaginary part of the lowest-order polarizations. Since the nuclear response 
is being probed in the spacelike (q 2 < 0) region, NN pairs can not be excited in these 
lowest-order descriptions. They can, however, be virtually excited and, thus, will become an 
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essential ingredient of the RPA response. Traditionally, NN excitations have been divided 
into two contributions, one being vacuum polarization and the other consisting of the Pauli 
blocking of NN excitations due to the filled Fermi sea [13| . The latter contribution is finite 
and has been shown to be essential for the conservation of the electromagnetic current. 
This contribution has been included in our calculations. The former, however, is divergent 
and must be renormalized. Since we are using a nonrenormalizable theory with derivative 
couplings, the renormalization of these divergent contributions becomes ambiguous at best. 
Thus, in order to avoid including ad-hoc parameters (e.g., cut-offs) we have decided to 
simply ignore the effect from vacuum polarization. Note that the (finite) real parts of the 
various polarizations have been calculated analytically and most of them have been published 
already ||23|| . Here we calculate them numerically so we can extend the formalism to include 
momentum-dependent self-energies. 



D. Momentum dependent vector and scalar self energy 

In a mean-field approximation to the Walecka model the vector (V) and scalar (S) self- 
energies are replaced by their classical expectation values. In this approximation the nucleon 
self-energy is real and energy independent. However, as the momentum of the nucleon 
becomes large there is an important coupling of the nucleon to nuclear excitations. Indeed, 
at intermediate energies it is known that the reactive content of the reaction is dominated 
by quasifree nucleon knockout. Thus, at large-enough momentum the nucleon self-energy 
will become complex and energy dependent. In order to calculate the nuclear response for 
a broad range of momentum transfers, we incorporate momentum-dependent self-energies 
into the nucleon propagator. Since a microscopic calculation of the energy dependence of 
the nucleon self-energy in the Walecka model awaits, we have used the phenomenological 
optical potentials of Ref. [[RJ . A detailed discussion of this momentum-dependent correction 
can be found in Ref. [ p~5|j . 

In a calculation with momentum-dependent self-energies, the effective mass and energy 
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of nucleons in the medium are no longer equal. In particular, the nucleon propagator must 
now be described in terms of Dirac spinors with masses and energies given by 



m; = m + s( P ) ; e p = e; + v( P ) = J P 2 + m; 2 + v( P ) 



(50) 



The basic formalism, however, remains unchanged except for the inclusion of a more realistic 
nucleon propagator given by 

1 



G(p) = {f + Ml) 



ITT 



, w> — + —5(p* -E*)6(k F -\p\) 



(51) 



where p* M = (p° — V(p),p) . Note that the vector potential can no longer be eliminated 
from the integrals defining the polarization by a simple change of variables. 



III. RESULTS 

In this section we present results for the inclusive cross section using a variety of approx- 
imations. We consider impulse-approximation calculations using, both, a relativistic Fermi 
gas of nucleons of mass M and a self-consistent ground state with an effective nucleon mass 
of M*. We also present two calculations including RPA correlations either with mass M or 
M*. The mass M RPA calculations can be directly compared to similar nonrelativistic cal- 
culations. Finally, results will be shown using momentum-dependent self-energies obtained 
from the phenomenological fit to the nucleon optical potential of Ref . 0] . We refer to this 
last calculation as "impulse with M*". The impulse with M calculation is commonly used 
to extract from experiment the mass parameter Ma present in the dipole fit to the axial 
form factor Ga- Our main goal is to estimate the sensitivity of this parameter to various 
relativistic nuclear-structure effects. 

The effective nucleon mass M* (in mean field theory) is obtained from a solution to the 



self-consistency equation [11 



a 2 A r k F M* 
M* = M-^-—^- d 3 k . (52) 
m 2 s (2vr) 3 Jo k 2 + M *2 
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Choosing the couplings to reproduce the bulk properties of nuclear matter at saturation 
[g 2 (M 2 /m 2 ) = 267.1] leads to a value of the effective nucleon mass of M* = 638 MeV at an 
assumed average density corresponding to kp = 225 MeV. 

In Fig. ^| we show the double differential cross section cPa j 'dVLk'dE^ [see Eq. ([12|)] for a 
neutrino energy of E v = 1 GeV. Fig. [3] (a) is for a momentum transfer of |q| = 0.5 GeV 
while Fig. ||| (b) is for |q| = 1.2 GeV. At |q| = 0.5 GeV the peak positions from the M* 
and M* calculations are shifted by less than 50 MeV relative to the Fermi-gas peak. This 
value represents an average binding-energy shift and has been observed experimentally in 
quasielastic electron scattering. Thus, a simple mean-field calculation is expected to give a 
reasonable description of the nuclear response. The situation changes considerably, however, 
at |q| = 1.2 GeV [Fig. |3| (b)]. Here, the M* calculation predicts a shift in the peak position 
that is substantially larger than the one obtained using the more realistic (M*) self-energies. 
Note that the impulse with M* calculation still predicts a substantial shift (of the order of 
50 MeV) relative to the Fermi-gas value. Also note the kinematical cutoff in q — above 
this value the scattering is prohibited kinematically. The impulse with M* calculation has 
a considerable amount of strength shifted into this inaccessible region. At this momentum 
transfer, the mean-field calculation overpredicts the binding energy shift seen in electron 
scattering [[0|. In contrast, the calculation using the impulse with M* shows a reasonable 



binding-energy shift and is practically insensitive to the kinematical cutoff. 

Fig. [| shows the inclusive RPA cross section at |q| = 0.5 GeV for various values of the 
Landau-Migdal parameter g'. For reference, the impulse results with M and with M* have 
also been included in Fig. |4] (c). The solid line shows the RPA response of a Fermi gas 
ground state while the dashed line shows the RPA response of the mean-field ground state. 
Since the pion contributes only through g', the softening and enhancement of the response 
in Fig. ^ (a) is exclusively due to an attractive — and thus unrealistic — rho-mediated 
residual interaction. As the value of g' is increased to g' = 0.3, the peak moves to higher 
excitation energy and the overall strength of the response is reduced. This is consistent 
with an additional repulsive component to the residual interaction arising from g'. Finally, 
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Fig. |] (c) shows the RPA results using the standard value of g' = 0.7. One now observes a 
slight quenching and hardening of the response due to the large value of g' . 

Since most of the experiments report the energy- integrated cross section da/dQ 2 we have 
taken our results for the double differential cross section and integrated over the allowed 
kinematical region of go 

da fQc it d 2 a 

dq , (53) 



dQ 2 Jo E u k> dE^dtt 
where the energy cutoff is given by 

q 2 — m 2 E v m 2 , 
AE V q A — m~ 

and m M is the mass of the produced lepton. In Fig. |5| (a) the energy-integrated cross section, 
da/dQ 2 , at E v = 1.2 GeV is shown for three calculations: impulse with M, RPA with M* 
using g' = 0.7, and impulse with M* Note that there is a substantial quenching due to RPA 
correlations at small momentum transfers, i.e., Q 2 < 0.3 GeV 2 . In the intermediate range 
(0.3 < Q 2 < 0.9 GeV 2 ) the importance of RPA correlations diminishes and no significant 
differences are observed between the three models. At an even larger Q 2 the RPA curve 
splits from the other two indicating the breakdown of the mean-field approximation [see 
Fig. |H (b)]. This breakdown, however, is not associated to RPA effects — which are no 
longer effective at this momentum transfer — but rather from the mean-field M* effects. 
Indeed, an RPA calculation using the free nucleon mass is within one percent of the impulse 
with M calculation at high Q 2 . 

We report similar calculations for antineutrinos in Fig. [| (b). The cross sections are 
considerably reduced relative to neutrinos because of the sign change in the vector-axial 
interference term [see Eq. (p3"|)]. In addition, since the cross sections fall to (almost) zero 
for Q 2 > 1 GeV 2 , most of the quasifree strength is located below the cutoff Q c . Hence, 
the kinematical cutoff does not have a significant effect on the M* curve at large Q 2 . At 
a smaller momentum transfer, Q 2 < 0.3 GeV 2 , and just as in the neutrino case, RPA 
correlations substantially reduce the cross section. At a larger Q 2 , most of the differences 
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observed between the RPA and the impulse with M calculations are due to the strong scalar 
potential. The impulse with M* calculation (dot-dashed curve) smoothly interpolates these 
two models. 

In the BNL experiment ||, the axial mass Ma, which controls the Q 2 falloff of the axial 
form factor [see Eq. ( |A9|) ], was extracted from antineutrino data using a Fermi gas (i.e., 
impulse with M) formalism. From an analysis of their data in the range Q 2 = 0.2 — 1 GeV 2 , 
the BNL group extracted a value of M A = 1.09 ± 0.03 (stat) ±0.02 (syst) GeV. In Fig. | 
we show the cross section, da/dQ 2 , averaged over the BNL neutrino spectrum. In Fig. ^ (a) 
we present four different calculations: impulse and RPA with M (solid lines), and impulse 
and RPA with M* (dashed lines). The two curves using the free nucleon mass M start to 
overlap at Q 2 = 0.4 GeV 2 and, insofar as we can regard them as our nonrelativistic limit, 
they agree with the nonrelativistic results obtained by Singh and Oset 0. Similarly, the 
two M* curves also coincide for Q 2 > 0.4 GeV 2 illustrating the fact that RPA correlations 
become unimportant in this Q 2 region. In Fig. |6] (b) we show three impulse curves in the 
Q 2 = 0.2 — 1 GeV 2 range — the region sampled in the BNL experiment. At Q 2 = 0.2 there 
is a 6 percent difference between the impulse with M calculation and the M* and M* values. 
This difference becomes larger, 20-30 percent, at Q 2 = 1.0 GeV 2 . It is, therefore, essential 
to estimate the sensitivity of Ma to these nuclear-structure effects. 

The raw experimental data as a function of the four-momentum transfer Q 2 was fitted 
with a theoretical curve to determine Ma- Since the experiment suffers from an uncertainty 
in the overall normalization, we use the ratio of the cross section at two different values of 
Q 2 to extract Ma- In Table |we show the ratio of the cross section using Q 2 = 0.2 GeV 2 and 
Q 2 = 1.0 GeV 2 . The value of the axial-mass parameter M A was varied in the 1.09 — 1.30 GeV 
range and the ratio of cross sections reported for various nuclear-structure models. For 
Ma = 1-09 GeV — the value extracted from the BNL data — the impulse with M ratio 
equals 19.42. This represents our baseline value for the ratio. In order to reproduce this 
value using the impulse with M* calculation the value of the axial-mass must be changed to 
Ma = 1-30 GeV (a 20 percent increase). Similarly, the value of Ma must be changed to 1.25 
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GeV in the RPA with M* and to 1.20 GeV in the M* calculation. There is, however, no 
Ma value within this range that can reproduce the ratio using an RPA with M calculation. 
This is because RPA effects substantially reduces the cross section at Q 2 = 0.2 (about 10 
percent) which was not expected from nonrelativistic calculations. 

Since we have established that RPA correlations become unimportant for Q 2 > 0.4 GeV 2 , 
we can eliminate the sensitivity to RPA effects by computing the ratio of cross sections using 
Q 2 = 0.5 GeV 2 and Q 2 = 1.0 GeV 2 (see Table 0). In this case, using the impulse with M 
calculation we obtain a baseline value for the ratio of 5.47 at Ma = 1-09 GeV. Now the 
RPA with M does not induce any change in the value of Ma — in agreement with the 
nonrelativistic calculation of Singh and Oset . This, however, is not the case for the other 
nuclear-structure calculations. Indeed, even in the impulse with M* calculation — which 
generates the smallest change — a 10 percent increase in the value of Ma is required. This 
10 percent uncertainty in the value of Ma could complicate the extraction of strange-quark 
matrix elements from experimental studies of neutral weak form factors [[/]] . 

IV. CONCLUSIONS AND OUTLOOK 

We have used a relativistic formalism to study inclusive charged-current neutrino scat- 
tering in an impulse approximation. The nuclear-structure information is contained in a 
large set of nuclear-response functions that were evaluated in nuclear matter using a variety 
of approximations. The simplest approximation that we considered was a relativistic free 
Fermi gas. This approximation was used by the BNL group to extract the axial mass param- 
eter Ma- We have used the Fermi-gas approximation — together with the BNL-extracted 
value of Ma = 1-09 GeV — to fix the Q 2 dependence of the inclusive cross section. We have 
examined the sensitivity of this extracted value for Ma to various nuclear-structure effects 
- such as the ones arising from the mean fields (M*), RPA correlations, and momentum- 
dependent self-energies. In essence, we have computed the changes in Ma that were required 
to reproduce the baseline free Fermi gas cross section once these additional nuclear-structure 
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effects were taken into consideration. This analysis is useful because a robust value of Ma 
is essential for a reliable determination of the strange-quark content of the nucleon. 

Our results indicate important corrections to the nuclear response due to the mean 
fields (M*) and to RPA correlations. Indeed, changes as large as 20 percent in Ma were 
observed whenever the whole range of Q 2 values (0.2 GeV 2 < Q 2 < 1.0 GeV 2 ) used in 
the BNL experiment were employed for the extraction. This uncertainty, however, can be 
substantially reduced. For example, it is well known that RPA correlations are effective 
only at small momentum transfers (i.e., Q 2 < 0.3 GeV 2 ). In addition, phenomenological 
fits to the nucleon optical potential indicate that, at large nucleon momenta, the mean 
fields are considerably smaller in magnitude than the ones predicted by the mean-field 
theory. Hence, some of the nuclear-structure uncertainties can be removed by employing 
phenomenological (momentum-dependent) mean fields and by restricting the range of Q 2 to 
the region 0.5 GeV 2 <Q 2 < 1.0 GeV 2 . Note that, although weaker than in the mean-field 
theory, the momentum-dependent optical potentials are still large and induce nontrivial 
changes in the effective mass and energy of a particle in the medium. Indeed, even in 
this best-case scenario a 10 percent uncertainty in M4 persists. With this extra 10 percent 
uncertainty, the BNL experiment |3J] by itself no longer provides strong evidence for a nonzero 
strangeness content of the nucleon 

In the future we will employ the present relativistic RPA formalism to calculate the 
inclusive cross section for atmospheric neutrinos. Atmospheric neutrinos have been observed 
over a wide range of energies (from a few hundred MeV to several GeV) in large water 
Cerenkov detectors — hence the need for a formalism that can address neutrino physics over 
this broad energy range. Atmospheric neutrinos are particularly interesting because of the 
current flavor anomaly in the z^j/V e ratio. This anomaly might signal neutrino oscillations. 
However, first one must rule out all conventional nuclear effects. Thus, it is important to 
examine these nuclear-structure corrections before a definitive statement about new physics 
can be made. 
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APPENDIX: 1 



We adopt the form factor parameterization used in Ref. M. First, the electromagnetic 
form factors are written in terms of simple dipole forms: 

G = (1 + 4.97t)~ 2 , r = -g 2 /4M 2 , (Al) 

Ff = [1 + t(1 + X p )]G/(1 + t), (A2) 

F 2 (p) = \ p G/(l + r), (A3) 

F[ n) = r\ n (l- V )G/(l + r), (A4) 

F 2 (n) = A n (l + rr / )G/(l + r). (A5) 
Here the anomalous moments are, 

\ p = 1.793, A n = -1.913, (A6) 

and 

rj = (1 + 5.6r) -1 . (A7) 

This parameterization is good for the neutron form factors provided r< 1. The isovector 
form factors are given by 

Fx = F[ p) - F} n \ F 2 = F 2 (p) - F 2 (n) . (A8) 



The axial form factor Ga is 



1.26 , 

" (1 - g2 /Mi)2 ' (A9) 



and the pseudoscalar form factor is given by 

2MG A 



, , , (A10) 
- q 2 
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APPENDIX: 2 



Analytic expressions for the imaginary parts of the polarizations introduced in Sec. [11 B 
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1,2,3) 



(A10) 



Ep = ^k 2 F + M* 2 
= Mm(Ep, E„ 
Emax = Max 



M*,E F -q ,-[-q + 1^1 



AM 



*2 



(All) 
(A12) 
(A13) 



The vacuum part does not contribute to the impulse response for spacelike momenta. 
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FIGURES 

FIG. 1. Feynman diagram for charged-current neutrino scattering. The neutrino with momen- 
tum k scatters off a nucleus in a state \ipi(pi)) and produces a muon or an electron with momentum 
k'. 

FIG. 2. A diagrammatic representation of the RPA correction AIIrpa. The RPA propagator 
Drpa is shown in (b) 

FIG. 3. Double differential cross sections for 1 GeV neutrinos. The solid curve is the momentum 
dependent calculation. Figure (a) is for |q| = 0.5 GeV while (b) is for |q| = 1.2 GeV. The dashed 
(dot-dashed) curve is the cross section obtained from a MFT (free Fermi gas) calculation. 

FIG. 4. Double differential cross section with RPA correlations for |q| = 500 MeV and E v = 1 
GeV. The solid curve is the RPA calculation with M while the dashed curve is the RPA with M*. 
(a) is for g'=0, (b) is for g'=0.3, and (c) is for g'=0.7. In (c) the M and M* impulse calculations 
are also shown for comparison. 

FIG. 5. da/dQ 2 for a neutrino (a) and antineutrino (b) scattering at an incoming energy of 1.2 
GeV. The solid curve is the impulse with M calculation while the dashed curve is the M* RPA with 
g' = 0.7. The effect of the momentum-dependent self-energies is indicated with the dot-dashed 
line. 

FIG. 6. da/dQ 2 averaged over the BNL spectrum. In (a) the two solid curves denote the 
impulse and the RPA (g' = 0.7) results using the free nucleon mass M. Also shown (two dashed 
lines) are the corresponding results using an effective nucleon mass M*. In (b) the cross sections 
are shown over the available experimental range in Q 2 . The solid and the dashed curves are the 
results of impulse calculations with M and M*, respectively. The dot-dashed is the calculation 
with momentum dependent self-energies. 
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TABLES 

TABLE I. The ratio of cross section ^p(Q 2 = 0.2 GeV 2 )/^(Q 2 = 1.0 GeV 2 ) for various 
nuclear structure effects at a given axial mass M A in the first column. 







^(Q 2 = 0.2GeV 2 )/^(Q 2 = 


1.0 GeV 2 ) 




M A (GeV) 


Impulse M 


Impulse M* 


RPA M* 


RPA M 


m; 


1.09 


19.42 


25.83 


23.72 


17.21 


21.85 


1.14 


18.47 


24.71 


22.59 


16.32 


20.83 


1.18 


17.75 


23.78 


21.69 


15.66 


20.03 


1.21 


16.91 


22.66 


20.61 


14.89 


19.08 


1.25 


15.95 


21.36 


19.38 


14.02 


18.00 


1.30 


14.9 


19.92 


18.04 


13.07 


16.81 



TABLE II. The same as Table | for the ratios -$ji{Q 2 = 0.5 GeV 2 )/^(Q 2 = 1.0 GeV 2 ). 







^(Q 2 = 0.5GeV 2 )/^(Q 2 = 


1.0 GeV 2 ) 




M A (GeV) 


Impulse M 


Impulse M* 


RPA M* 


RPA M 


m; 


1.09 


5.47 


6.50 


6.40 


5.34 


5.77 


1.14 


5.35 


6.38 


6.27 


5.21 


5.65 


1.18 


5.24 


6.26 


6.14 


5.10 


5.54 


1.21 


5.10 


6.09 


5.96 


4.96 


5.39 


1.25 


4.92 


5.88 


5.75 


4.78 


5.21 


1.30 


4.72 


5.62 


5.49 


4.58 


4.99 
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